PREDICTION OF CAVITY SIZE IN THE PACKED BED SYSTEMS 
USING NEW CORRELATIONS AND MATHEMATICAL MODEL 

Field of the Invention: 

The present invention relates to prediction of cavity size in the packed bed systems using 
new correlations and mathematical model. Simplified equations, based on analytical 
solution of one-dimensional mathematical model, have been developed along with the 
cavity correlations to describe the cavity size and hysteresis. The proposed correlations 
and mathematical model give a universal approach to predict the cavity size which is 
applicable to any packed bed systems like blast furnaces, cupola, Corex, catalytic 
regenerator, etc. and is able to represent, in a good way, the data of other researchers 
provided the frictional properties of the particulate are known. Developed correlations 
and model can be used directly to optimize the above mentioned and other related 
processes. 

Prior Art: 

On Packed Bed: In the packed bed, contact forces between the particles and wall-particle 
have been considered widely in explaining its behavior in various conditions. Reference 
may be made to F.J. Doyle III, R. Jackson and J.C. Ginestra, "The phenomenon of 
pinning in an annular moving bed reactor with crossflow of gas", Chem. Eng Sci, 41(6) 
1986 1485, wherein they have studied moving bed of cross flow theoretically in order to 
study the pinning effect in catalytic reformer. Their analysis is based on force balance 
approach considering the gas drag, stresses and gravity forces. The drawbacks of their 
simplified model, which they had presented, are (i) it is based on arbitrary assumption of 
the radial variation of the stress in the moving beds. Due to this reason their numerical 
values are greater than limited experimental values by a factor of two. (ii) They had 
assumed the shear stress at the wall of the moving bed reactor to act in the downward 
direction, (iii) The analysis was confined to the growth of cavity till the solid flow ceases 
in moving bed. 

Reference may be made to V.B. Apte, T.F. Wall and J.S. Truelove: AIChEJ, 1990, vol. 
36 (3), pp. 461-468, wherein they have analysed the stress distribution above a cavity 
formed by an upward gas blast from the bottom of a two-dimensional packed bed. They 
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wrote one dimensional elemental force balance, along the streamline coincident with the 
tuyere axis, between the pressure, bed weight and frictional forces. The drawbacks o their 
model are (i) they had assumed that frictional stresses always act in the upward direction, 
(ii) They were unable to show any hysteresis results, (iii) They neglected any acceleration 
effect due to slowing down of the gas and (iv) did not predict the cavity size. Mainly their 
study was concentrated on the stress distribution in the packed bed under increasing 
velocity. 

Reference may be made to J. F. MacDonald, and J. Bridgwater, Chem. Eng. Sci., 1997, 
vol. 52 (5), pp. 677-691, wherein they have studied the phenomenon of void formation in 
stationary and moving beds of solids and unified the behaviour using dimensional 
analysis. The drawback of their correlation is that they recognised the importance of 
frictional forces in cross flow but were unable to include it in their dimensional analysis. 

On (Ironmaking, Lead, Corex, etc.) Blast Furnaces: In the blast furnace, gas is introduced 
laterally at a high velocity through a pipe, called tuyere, in the packed bed of coke. This 
creates a cavity in front of the tuyere called raceway. Coke is burnt in this zone to supply 
heat to the process. Therefore, coke particles get consumed in this region and they are 
replenished by fresh coke particles from the top of the raceway. So the whole burden 
descends in the downward direction. The size and shape of the raceway affects the 
aerodynamics of the furnace and thus affects the overall heat and mass transfer. Due to 
this reason, raceway has been studied extensively both theoretically and experimentally. 
In case of blast furnace, many authors have presented raceway correlations to predict the 
raceway size which are listed in Table 1. Most of these correlations are based on cold 
model study and some of them are based on hot model and plant data study. 

References may be made to J.D. Lister, G.S. Gupta, V.R. Rudolph and E.T. White: 
CHEMECA'91 Conf., 1991 Newcastle, Australia, vol. 1, 476 and S. Sarkar, G.S. Gupta, 
J.D. Litster, V. Rudolph, E.T. White and S.K. Choudhary: Metall Trans,, 2003, 34B (2), 
183-191, wherein they have that none of these correlations predicts the raceway size in 
industrial conditions reasonably and they also differ to each other. It is observed that all 
the experimental correlations have been based on various forms of Froude number. The 
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raceway size has been institutively correlated with this number along with some other 
parameters such as height of the bed, width of the model and tuyere opening. 

References may be made to J.F. Elliott, R.A. Bachanan and J.B. Wagstaff: Trans. AIME, 
1952, vol. 194, pp; 709-717. J. Taylor, G. Lonie and R. Hay: JISI, 1957, vol. 187, p330; 
J.B. Wagstaff and W.H. Holman: Trans. AIME, March 1957, pp. 370-376. M; Hatano, B. 
Hiraoka, M. Fukuda and T. Masuike: Int. ISIJ, 1977, 17, pp. 102-109; M. Nakamura, T. 
Sugiyama, T. Uno, Y. Hara and S. Kondo: Tetsu-to-Hagane, 1977, vol 63, pp. 28, wherein 
one can see that these correlations (see Table I) are not evolved based on a systematic study 
i.e. by applying dimensional analysis and finding the relevant groups. 

On the other hand, theoretical correlations have been obtained by simplifying the actual 
theoretical equations logically by P.J. Flint and J.M. Burgess: Metall. Trans., 1992, vol. 
23B, pp. 267-283 and J. Szekely and J.J. Poveromo: Metall. Trans., 1975, vol. 6B, pp. 119- 
130. These correlations are more systematic. Also, all the empirical correlations, for the 
two and three-dimensional models, have been obtained for the velocity increasing case. 

It must be mentioned here that one can get two raceways size at the same gas velocity 
depending on whether the measurement is made in the increasing or decreasing gas 
velocity. This phenomena is called raceway hysteresis. References may be made to J.D. 
Lister et al. 1991 and S. Sarkar et al. 2003, wherein hysteresis phenomenoa has been 
described in detail and has been reported that the decreasing velocity correlation is more 
relevant to blast furnace. 

Since the raceway size in the increasing and decreasing velocity case vary by 
approximately a factor of 4, the raceway size can affect considerably the predictions of 
heat, mass and momentum transfer in the blast furnace. At this juncture something about 
the raceway hysteresis should be mentioned because the background of the 
correlations/mathematical model developed in this study is based upon this phenomena. 
Reference may be made to S. Sarkar et al. 2003, wherein they have explained raceway 
hysteresis phenomenon in details and have proposed that raceway hysteresis can be 
represented by the following equation, based on their experimental results. 
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Pressure Force - Bed Weight ± Frictional Forces (Stresses) = 0 (1) 

The physical interpretation of this equation is that when the raceway is expanding, the 
particles near and above the raceway are being pushed in the upward direction. So the 
frictional stresses will tend to oppose this motion of the particles and hence act in the 
downward direction and is fully mobilized. When we start to decrease the blast velocity 
from a maximum value, the particles above the raceway are trying to fall down. So the 
frictional forces act against this movement and start increasing in magnitude in the 
upward direction progressively. Once the frictional stresses acting in the upward direction 
become fully mobilized, further reduction in blast velocity results in decrease in the 
raceway penetration. A positive sign in the equation (1) in the frictional forces term 
indicates cavity wall friction acting upwards (for velocity decreasing) and a negative sign 
indicates cavity wall friction acting downwards (for velocity increasing). Pressure force 
always acts in the upward direction and bed weight always acts in the downward 
direction. 

Objects of the Present Invention 

The main object of the present invention is to provide a method and a system for prediction 
of cavity size in the packed beds using new correlations and / or mathematical model 
which obviates the drawbacks as detailed above. 

Statement Of Invention : 

Accordingly the present invention provides a method and a system for prediction of 
cavity size in the packed beds using new correlations and mathematical model which 
comprises the development of two correlations, one each for increasing and decreasing 
gas velocity respectively based on n - theorem for two-dimensional cold model 
experiments having the variables like bed height, tuyere opening, void fraction, frictional 
and physical properties of various materials, gas flow rates and width of the model as 
well as it comprises the development of one dimensional mathematical model based upon 
a force balance approach (as discussed in prior art) and then solving the developed 
equations analytically for pressure force, frictional force and bed weight to describe the 
cavity hysteresis and to predict the cavity/raceway size & minimum spouting velocity/ 
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instability in packed beds and later on to compare the correlations and model results with 
experiments and published /plant data on cavity size. 

In an embodiment of the present invention it clarifies the direction of frictional forces and 
gives a logical explanation of it to describe the hysteresis in the packed beds. 

In another embodiment of the present invention it also brings out that decreasing velocity 
data are relevant to operating blast furnaces. 

In yet another embodiment of the present invention it gives, through math, model, the 
maximum operating gas velocity in a packed bed beyond which it will become unstable. 

Detailed Description of the Invention: 

Accordingly, the present invention provides a computer based method for determining 
the cavity size in packed bed systems using correlation or mathematical model, said 
method comprising the steps of: 

(a) obtaining data related to material properties of the packed bed system; 

(b) calculating the cavity radius for both increasing gas velocity and decreasing gas 
velocity using mathematical model incorporating the stresses/frictional forces as: 
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respectively; or calculating the cavity radius for both increasing gas velocity and 
decreasing gas velocity using mathematical equations based on correlation as: 
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respectively, and 

(c) calculating the cavity size using the cavity radius obtained in step (b). 

In an embodiment of the present invention, the data related to material properties of the 
packed bed comprise bed height, tuyere opening, void fraction, wall-particle friction 
coefficient, inter-particle frictional coefficient, gas velocity, model width and particle 
shape factor. 

In another embodiment of the present invention, the data related to the material properties 
of the packed bed include experimental data already obtained or on-line data. 

In yet another embodiment of the present invention, the frictional force (F wd ) in equations 
28 and 29 is given by: 
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In still another embodiment of the present invention, wherein to determine the cavity 
radius using increasing velocity correlation as given by equation 33 was developed using 
7r-theorem to get the important dimensionless numbers 
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where, symbols are Blast furnace radius W, Effective bed height H, Blast velocity v b > 
Tuyere opening D t , Void fraction 8, Gas viscosity |i g , Particle size d p , Shape factor (|) s , 
Density of gas p g , Density of solid p s , Coefficient of wall friction |j, w , acceleration due to 
gravity g, the effective diameter of the particle is given by d eff = d p <fi s , effective density 

of the bed is given by p eff = ep g + (l — s)p s , wall-particle frictional coefficient is given 
by ju w = tan (f) w , where, (|) w is an angle of friction between the wall and particle D r is 
cavity diameter and all units are in SI. 

In one more embodiment of the present invention, wherein to determine the cavity radius 
using decreasing velocity correlation as given by equation 36 was developed using 7i- 
theorem to get the important dimensionless numbers 
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where, symbols are Blast furnace radius W, Effective bed height H, Blast velocity Vb, 
Tuyere opening D t , Void fraction s, Gas viscosity (i g , Particle size d p , Shape factor (|> s , 
Density of gas p g , Density of solid p s? Coefficient of wall friction |i W5 Acceleration due to 
gravity g, the effective diameter of the particle is given by d eff =d p <f> s , effective density 

of the bed is given by p eff - ep g + (l — s)p s , wall-particle frictional coefficient is given 
by /u w = tan^ w , where, (|) w is an angle of friction between the wall and particle D r is 
cavity diameter and all units are in SI. 

In one another embodiment of the present invention, wherein the packed bed systems 
include blast furnaces, cupola, corex, catalytic regenerator. 
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It is important in any gas-solid process to achieve a uniform gas and solid distribution 
which determines its performance. Packed, spouted and fluidized beds fall under these 
categories and are widely used in industries. A common feature of all these beds is that 
they all show hysteresis. Figure 1 shows a cavity hysteresis plot between the cavity 
diameter and gas velocity which clearly shows the presence of hysteresis. It is evident 
from the figure that cavity size increased with increasing gas velocity. When the gas 
velocity was decreased from the maximum value (A), there was initially almost no 
change in the cavity size. However, when a critical velocity was reached (B), the cavity 
size began decreasing with decreasing velocity, but was always larger than that for the 
same velocity achieved in increasing velocity. This is the cavity hysteresis phenomenon. 
The cavity hysteresis found in the packed beds is similar to hysteresis found in 
fluidization beds. 

Here we are presenting a one-dimensional theoretical model based on equation (1) to 
predict the cavity size and to describe the mechanism of hysteresis in the packed bed. 
Also we are presenting new cavity/raceway size correlations using 71-theorem. The 
various terms in the equation (1) are expressed in their mathematical form below. 

Model Formulation 

Let us consider a two-dimensional packed bed of solids of height H and width W as 
shown in figure 2. The gas was injected laterally at a particular blast velocity v b through a 
slot type nozzle of opening D T , creating a void of equivalent radius R in front of it. Let p 
and \x be the density and viscosity of gas respectively. d p is the particle diameter and 8 is 
the void fraction of the bed. D. Akamatsu, M. Hatano and M. Takeuchi, Tetsu-to-Hagane 
58 (1972) 20, had measured the pressure inside the cavity and had found that that the 
pressure distribution was relatively uniform. Therefore, it is reasonable to assume that the 
gas flows radially from the center of the cavity into the surrounding packed bed with the 
velocity varying along concentric circles. Based on several other experimental and 
theoretical studied (Szekely & Poveromo (1975), Flint & Burgess (1992), V.B. Apte, T.F. 
Wall and J.S. Truelove, Gas Flows in Cavities Formed by High Velocity Jets in a Two- 
Dimensional Packed Bed, Chem Eng Res Des 66 (1988) 357, and, M. Hatano, K. Kurita 
and T. Tanaka, Ironmaking Proc. Iron Steel Soc 42 (1983) 577), isobaric condition inside 
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the cavity has been assumed. The velocity of gas, which is moving upwards, varies with 
the r-direction (distance from the center of cavity) but does not vary in the angular 
direction. 

Pressure exerted by the gas: It has been reported (Flint & Burgess, 1992 and Apte et. al., 
1990) that gas velocity becomes almost constant at the exit bed velocity after some 
distance say r = r D from the cavity center. The corresponding velocity at this distance is v 
= v H (see fig. 2). Then, equating the mass flow rate of gas at the nozzle and at a distance 
r Q from the center of the cavity, we get 

p v b D T = p (27ir 0 - D T ) v H Or, v H = v b D T / (27ir 0 - D T ) (2) 
Also, equating the mass flow rate of blast at the nozzle and at the bed surface, one gets 

p v b D T = p W v H Or, v H = Vb D T / W (3) 
From (2) and (3), one obtains 

r G = (W + D t )/2ti (4) 
After a distance r Q from the cavity center, the velocity of the gas will be constant. This 
observation has been confirmed computationally by Flint & Burgess, 1992. Analysis of 
the experimental data of Apte et al., 1990, also verifies the validity of equations (3) & (4). 

Let v(r) be the gas velocity at a distance r from the center of the cavity. Then on equating 

the mass flow rate at the nozzle opening and at a distance r from the center of the cavity, 

p (27ir - D T ) v(r) = pv b D T Or, v(r) = v b D t /(2tit - D T ) (5) 

Based on the above equations, the modeled velocity profile may be written as 

v(r) = VbD T /(27tr-D T ) ,r<r 0 

= v H , r > r Q (6) 

Therefore, the velocity of the gas, which is moving upwards, varies inversely with 
distance up to a distance of r G from the center of the cavity (radial region) and then 
remains constant beyond this (cartesian region). 

For drag force in fluidized bed, many researchers have widely used Richardson-Zaki 
correlation. Similarly, in packed bed the force per unit volume exerted by the gas on solid 
is given from well-known Ergun's equation 
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-^- = a V (r)+/3v 2 (r) (7) 
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where ' a = sw; "* A= 

<|) s is the shape factor of particle. In practice, the gas velocity in the radial region is high. 
At these high velocities the viscous term is negligible compared to the inertial term i.e. 
av(r) « p\ 2 (r) . Therefore, the force exerted by the gas on the solids is given by 
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Similarly, the force exerted by the gas in the cartesian region would be 

F 2 = \-&- ( W + D A dz where, -dpi dz = av H + pv 2 H (9) 
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And, (W+D T )/7i = (2r Q ) is the diameter of the largest circle, in the varying velocity 
region, through which the gas flows out radially and enters into the cartesian region as 
shown in Figure 4. z is the variable height of the packed bed from the tuyere level. After 
integrating the equation (9), one gets 

F 2 = (cc+p v H ) v h [(W+D t )/tc] [H- (W+D t )/2ti] = (cc+p v H ) v H (2r Q ) (H-r 0 ) (10) 

Therefore, the total force exerted by the gas (either in increasing or decreasing velocity) 
on the solids above cavity can be given by 

F pr - f =F,+F 2 

Determination of Frictional Force in the Cartesian Region (Decreasing Velocity): In 

decreasing velocity, the particle-wall frictional force acts in the upward direction as 
explained earlier and is shown in figure 3 along with other forces in which z-axis is along 
the upward direction from the tuyere level or from the center of the cavity. It is assumed 
that the normal stress (cr 2 ), acting in the upward direction, is constant at any distance z 
from the bed surface, dz is the thickness of the slice over which the elemental balance is 
done. Let a 2 + da z be the reaction stress at distance z + dz acting in the downward 
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direction and z w be the particle-wall frictional stress. M is the bed weight per unit 
volume. Equating the forces acting on the element, we get: 

(a z + da z )xWx\ + MxWxdzxX = a z xWx\ + 2r w xdzx\ + dPxWxX (11) 

Factor 2 in the second term on the right side is due to r w acting on both sides of the wall. 
dP is the force per unit area exerted by gas over the element = (-dp/dz) dz. 

Following Janssen approach (H.A. Janssen, Versuche uber getreidedruck in solozellen. 
Ver. Deutsch. Ing. Zeit. 39 (1895) 1045), it is assumed that the vertical stress (cr z ) and 
horizontal stress (cr x ) are the principal stresses. Therefore, particle-wall frictional stress 
can be written as t w = |^ w Ka z . Where, K = ((l-sin(J>)/(l+sin(|))) is the lateral pressure 
coefficient K and § is the angle of internal friction. |j, w is the coefficient of friction 
between the bed walls and the particle. Substituting the value of t w in the equation (11) 
and after some simplification one gets 

d^_2M^ = _ M _dp 

The solution of equation (12), using the boundary condition, at z = H, cr z =0, would be 

a z =f{l-e^}-e-)^e^d Z (13) 
C ^ dz 

where, C = 2|a w K/W, is the bed support factor. The first term on the right hand side of 
equation (13) is the effective bed weight, while the second term represents the upward 
gas pressure drag. For a uniform gas flow in the bed i.e. a constant -dp/dz (= a v H + P 
v H 2 ), equation (13) reduces to (after substituting the value of v H from the equation (3)) 

For a no gas flow situation, i.e. a static bed, the equation (14) reduces to 

{l-e C( -">} (15) 
This is a classical Jansen's equation, assuming a constant cr z over any horizontal cross 
section. For deep beds as (H-z) oo , the above equation becomes a z = M/C. 
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C is a function of W, // w and K and hence is a measure of the particle-wall frictional 

support. Larger C implies a larger particle-wall frictional support and hence a smaller 
effective bed weight. Also C is inversely proportional to the width of the model. Larger 
the width of the model, lower will be the value of C. From equation (15), 
as lim c _> 0 cr z = M(H - z) , implying that for C=0, the bed weight would be transmitted as 
an equivalent hydrostatic head. 

It is necessary to determine the particle-wall frictional force acting over this region in the 
cartesian system. The particle-wall frictional force F wd2 acting in the upward direction in 
the region lying over a distance 2r G , in the constant velocity region is obtained by 
multiplying particle-wall frictional stress t w by area and integrating it from z= r 0 to z = H. 

K«2 = K><2 1^ - T J dz = \/u w Ka z x 2 [ ^ j dz 



= W 



Frictional Forces in the Radial Region: Like the cartesian region, the radial system for 
elemental balance is shown in figure 4. Resolving all the forces along the radial direction 
and doing a force balance over the top portion of the circular element, one gets 

<j r {n{27rr-D T )} + 2r w xdrxl + dPx {n(27tr - D T ) }xl 
= (cr r + da r ){n(2x(r + dr)-D T )} + hx M {n{2nr - D T ) }xdrxl 

where, dr is the thickness of the circular section over which the elemental balance is 
carried out. a r is the radial stress at radius r and cr r + da r is the reaction stress at radius 
r+dr. x w is the particle wall frictional stress acting in the upward direction, n is the factor 
of contribution of top portion of the cavity to the total cavity area and h is the factor 

arising due to resolving the vertical force along the radial direction, so h = Jcosor da 

dP=-^xdr and -f = fr'fr) = , ^ , 2 (18) 
dr dr (2nr-D T ) 
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Assuming that a r and a Q are the principal stresses, then r w = ju w cr e = ju w Ka r (19) 
After substituting the value of x w and dP in the equation (17) and integrating it, one gets 
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Close to the cavity region, where the velocity is very high, the second term in the above 
equation becomes significantly high leading to a drop in stress. A, the constant of 
integration, can be calculated using the boundary condition at the interface of the radial 
and cartesian systems i.e. at r = r c , cr r = cr z . Finally, equation (20) can be written as 
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In the above equation, the bed weight (M) containing terms when added give the 
effective bed weight and the blast velocity (v b ) containing terms when added give the 
effective upward gas pressure drag. The wall frictional force can be obtained by 
multiplying t Wj resolved along the vertical direction, with the area and integrating from r 
= R to r = r Q as given below. 

F**x = 2p jV w {n(2nr) - D T }dr = 2p )» w K(j r {n(2nr) - D T }dr (22) 

R R 

where, p-h = factor obtained by resolving the radial force along vertically upward 
direction. On integration, equation (22) can be written as 
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Elemental Force Balance During Increasing Velocity: 

It can be done in a similar way as it has been done for the decreasing velocity case as 
reported in S. Rajneesh, M.E. (Int.) Thesis, Indian Institute of Science, Bangalore, Sept. 
2000 and in CSIR Report No. 22(285)/99/EMR-II. 

Force balance over top of the cavity (for Decreasing Velocity) 

From equation (8), force exerted by the gas above the top of the cavity in varying 
velocity region in the vertically upward direction (after resolving it) is given by 
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Therefore, the total force exerted by gas on solid in the upward direction is 

F p , f =F ia + F 2 = jln|l-ln(*-|?-)j + (a + p\ H ) (H - r Q ) v H (lr 0 ) (25) 

Similarly, the total particle-wall frictional force acting in the upward direction is 



F w d — F\vtd\ + ^wd 2 



(26) 



where, Fwtdi and F w d2 are given by the equations (23) and (16) respectively. 



It is assumed that bed weight is transmitted hydrostatically over the cavity roof. For 
simplification it has been assumed that the contribution of the bed weight from the sides 
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to the cavity formation is negligible. Therefore, bed weight at the top of the cavity roof is 
= Bed weight/area x Area of the top portion the cavity = M{H - R) x n{2nK) x 1 (27) 

Wherein "n" is the factor of contribution of the top portion of the cavity to the cavity 
area. 

After substituting all forces (equations 25, 26 and 27) in the equation (1) and after some 
simplification one can write the equation in terms of cavity radius, R. 

(28) 

Solving equation (28) for R numerically, gives the cavity radius in the velocity 
decreasing case and thus the cavity diameter D r = 2R. 

Similarly, one can establish the force balance over the cavity in the case of increasing 
velocity and can obtained the cavity diameter as explained above. 

(29) 

Raceway/Cavity Size Correlations: 

Velocity Increasing Case (Using BUCKINGHAM n - THEOREM) 
The raceway is formed due to a balance between the pressure force exerted by the gas, 
bed weight and the frictional forces as described by the force balance equation (1). The 
pressure force exerted by the gas comprises the inertial and viscous force. The inertial 
force exerted by the gas depends on the blast velocity (v b , m/s), density of the gas ( p g , 

kg/ m 3 ) and the tuyere opening (D T , m). The viscous force exerted by the gas depends on 
the viscosity (ju, Pa.s) of the gas and the particle diameter (d p , m). The bed weight 

exerted by the packing depends on the density of the solid ( p s , kg/ m 3 ), acceleration due 

to gravity (g, m/sec 2 ), height of the bed (H, m) and void fraction of the bed. The 
frictional forces (or stresses) depend on the internal and wall angle of friction and this 
causes the introduction of the wall-particle frictional coefficient ju w and v, the inter- 
particle frictional coefficient. Finally, the width of the bed W has been taken since it has 



15 



been varied during the experiments as it affects the raceway penetration. In other words, 
the raceway diameter (D r , m) in a packed bed is a function of the property of material 
used for packing, property of the gas injected through the tuyere, the geometrical 
parameters and the frictional parameters i.e. 

The effective diameter of the particle is given by d eff =d p sh, where d p = diameter of 

the particle and sh = shape factor of the particle. Effective density of the bed is given by 
p eff = sp g + (l - s)p s . Wall-particle frictional coefficient is given by ju w = tan <f> w and 

inter-particle frictional coefficient is given by v= tan^ . Where, (j> and (j) w are an 

internal angle of friction between the particles and angle of friction between the wall and 
particle respectively. 

D r = f (PeffrPg> V b>g> d eff>M>E>T> H > W >M W >V) ( 30 ) 

Since the total number of variables is 12 and the number of independent variables in 
terms of which the variables can be expressed is 3, the number of dimensionless groups 
that will be obtained from Buckingham n— theorem is 9. Using n- theorem, the 
correlation for the raceway diameter was obtained as: 

A group involving d p and D T has been omitted as some other groups already represent 
these quantities. Similarly, v has been neglected because in 2D cold model wall particle 
friction would be dominating rather than inter-particle friction. Moreover, the value of § 
changes with the gas flowrate (R. Jackson and M.R. Judd, Further consideration on the 
effect of aeration on the flowability of powders, Trans IchemE, 59 (1981) 119) which 
makes difficult to assign it single value. 

The first dimensionless group on the right side is related to pressure drop. Second group 
is Froude number which gives the ratio of inertial to gravitational forces. It is used to 
describe the gas/solid/liquid systems. Many previous authors have correlated raceway 
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size with this number. The third group is well known Reynolds number. The left hand 
side group of equation (30) is known as raceway penetration factor. 



From the experimental values, obtained in the velocity increasing case, the values of 
dimensionless groups are evaluated. The resulting data is then subjected to regression 
analysis to determine the constants a, b, c, d, e, f, and k. The values of the constants 
obtained are a=0.79, b=0.81, c=0.0035, d=0.88, e=0.89, £=-0.24 and k=243.5. 



From these values it is clear that Reynolds number is of least significance. All other 
parameters are important. Therefore, after neglecting the Reynolds number term and 
performing regression analysis again one gets the values of the coefficients as a=0.79, 
b=0.81, d=0.85, e=0.88, f==-0.23 and k=247. It can be observed that there is not much 
change in the value of the coefficients after neglecting the Reynolds number. The effect 
of the Reynolds number is negligible because of the inertial conditions prevailing during 
the raceway experiments performed. Since the values of coefficients a, b, d and e are 
quite close, we can group them into a single dimensionless group and the simplify form 
of the correlation can be written as: 



= k\ 



2 r\2 

T 



D T {p eff gd eff HW 



(32) 



Doing regression analysis again, we get the values of the coefficients as a=0.80, b= -0.25 
and k=164. The R 2 value of the correlation was found to be 0.96. Therefore, the final 
form of the correlation for increasing velocity is: 



= 164 



/ o t \0.80 

' pA d t 



yp eJf gd eff Hw^ 



CO 



-0.25 



(33) 



Velocity Decreasing Case: The correlation for the raceway diameter as before is given 
by: 



Pg 



D T ^Peff 



,2 W 



gdeff 



PgVbdeff 



c ' Dt W 



w 



(34) 



A regression analysis was performed on the experimental data, obtained in the gas 
velocity decreasing case, to determine the constants a, b, c, d, e, f, and k. The values of 
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constants obtained are: a=0.60, b=0.62, c=-0.024, d=0.51, e=-0.095, £=-0.235 and 
k=3.3612. The R 2 value of the correlation was found to be 0.96. 



As before, one can neglect the Reynolds number since its coefficient c is very small. 
Since the values of other coefficients a, b, and d are quite close, one can group these 
dimensionless groups into single group. Thus the simplify form of the correlation can be 
expressed as: 



where k, a, b and c have to be determined by regression analysis again. Using the above 
equation and performing regression analysis, one obtains the following final form of the 
correlation for decreasing velocity. 



Equations (32) and (35) are the desired raceway size correlations for increasing and 
decreasing velocity respectively. It is interesting to note that bed height and tuyere 
opening play an important role in increasing than decreasing velocity. The results 
obtained from these correlations will be compared with the experiments and plant data. 

Experimental Plan: 

Before the experimental procedure is described, it is necessary to distinguish between 
two types of two-dimensional apparatus (G.S.S.R.K. Sastry, G.S. Gupta and A.K. Lahiri, 
Ironmkg & Steelmkg, 30 (1) (2003) 61) which have been used by various researchers. 
These are classified as pseudo two-dimensional and two-dimensional models. In two- 
dimensional model a tuyere in the form of a rectangular slot is introduced across the 
entire thickness of the model. This ensures a uniform blast velocity across the entire 
width and no expansion of the jet in the third dimension takes place. Thus the 
phenomenon is confined strictly to two dimensions. In pseudo two-dimensional models, 
jet of air is introduced through a tuyere (mostly circular) placed in the longitudinal central 
plane of the model and the phenomenon is observed from the sidewalls where the effects 




(35) 




(36) 



The R 2 value of the correlation was found to be 0.96. 
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are visible. The jet can expand in front of the tuyere in all directions but it is assumed that 
there is negligible effect due to the jet expansion in the direction perpendicular to the 
tuyere axis. Most of the investigations on raceway have been done on pseudo two- 
dimensional model except by Flint & Burgess (1992), Litster et al. (1991) Sarkar et al. 
(1993), and (G.S.S.R.K. Sastry, G.S. Gupta and A.K. Lahiri, Int. ISIJ, 43 (2) (2003) 153). 
It is obvious that only two dimensional model can give better accuracy, therefore, only 
two dimensional models have been used in the present study. 

As such the raceway size is a function of physical and frictional properties of the material 
and geometrical parameters of the experimental setup. Therefore, many experiments were 
performed to obtain the raceway size as a function of these parameters in both increasing 
and decreasing gas velocity. Table 2 shows the range of various variables (geometrical) 
along with experimental variables used during the experiments. All the particles, which 
were used during the experiments, were having the ratio of apparatus thickness (opening) 
to particle diameter always greater than 12 or more in order to avoid the wall effect. All 
experiments were carried out in two-dimensional cold models which were reinforced 
using iron bars to prevent the bulging. PVC slot tuyeres were used. A schematic diagram 
of the equipment is shown in Figure 5. 

The bed was packed with a desired material to a desired bed height above the tuyere 
level. Room temperature air was used as the blast gas to form the raceway. The air flow 
rate to the tuyere was increased gradually until the point at which the raceway just began 
to form, then it was shut off immediately. This procedure was necessary to clear the 
tuyere of the beads which entered the tuyere when the bed was filled. The air flow rate 
was then increased gradually from zero to the fluidisation limit of the bed in steps. At 
each step, two minutes were allowed for the raceway size to reach equilibrium, then the 
raceway penetration (size in the gas entry direction) and height were measured directly 
using a ruler and tracing the raceway boundary on a transparent graph paper. When the 
maximum gas flow rate for the experiment was reached, the flow rate was reduced 
through the same steps. Raceway penetration and height were measured in thesame way. 
Each experiment was repeated at least thrice. However, average value has been reported. 
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Various physical properties of the materials used in the experiment, are listed in Table 3. 
Hundreds of experiments were performed to obtain the raceway size by changing the 
dimensions of the apparatus, bed height, tuyere opening, gas flow rate and material 
properties. 

Table 2: List of geometrical and experimental variables 



Apparatus 
Number 


Bed dimensions 
(HXWX T), 
mm 


Tuyere 

opening 

(mm) 


Gas 

velocity 
(m/s) 


Bed 

height, 

m 


Material 


Experimental 
condition 


1 


2300x1000x100 


6,10,25, 
50&79 


0-120 


0.2-1 


Polyethylene 


Both 

(increasing and 

decreasing 

velocity) 


2 


1800x600x60 


5 


0-110 


0.2-1 


Glass, Plastic 


Both 


3 


830x380x40 


5.5 


0-40 


0.1-0.5 


Plastic, 
Mustard seed 


Both 


4 


700x285x17 


5 


0-25 


0.1-0.5 


Quartz 


Increasing 



Table 3. 



Physical properties of the materials 



Material 


Shape 


Density 
(kg/m 3 ) 


Particle 

diameter 

(mm) 


Particle 

wall 

friction 

(Hw) 


Shape 
factor 


Min. 

fluidization 
velocity (m/s) 


Void 
fracti 
on 


Plastic 


Spherical 


1080 + 20 


5.8 + 0.04, 
2.1 ± 0.1 


0.22 


1.0 


1.37 (5.8mm) 
0.67 (2.1mm) 


0.42 


Polyethy- 
lene 


Cylindrical 


920 ± 30 


4.1 (Equiv. 
Dia.) 


0.29 


0.87 


0.84 


0.42 


Glass 


Spherical 


2770± 90 


2.7± 0.01 


0.16 


1.0 


1.39 


0.43 


Quartz 


Irregular 


2550 ±70 


Equiv. Dia. 
1.09,1.55mm 


0.2 


0.65 


0.87 (for 
1.55mm) 


0.4 


Mustard 


Spherical 


1070 ± 10 


2.2 ± 0.2 


0.22 


1.0 


0.69 


0.39 



Scientific Explanations: 

Below, numerical results have been presented, based on developed mathematical model 
here, considering an experimental apparatus number 1 and polyethylene beads (see 
Tables 2 & 3). The angle of wall friction and inter-particle friction were measured using 
shear apparatus which were 15.6 and 38 respectively. However, in order to know the 
angle of inter-particle friction in presence of air, the equation suggested by Jackson and 
Judd (1981) was used which gives the value of internal single of friction 28 at an average 
gas velocity of 40m/s. Height (H) of the packed bed above the tuyere level is lm. The 
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total width (W) and thickness of the model are lm and 0.1m respectively. The value of r 0 
[=(W+D t )/27c] is 0.16m. Therefore, the system is cartesian from 0.16m to lm (top of the 
bed surface). Value of n, measured experimentally, was 0.8. Before comparing 
experimental, plant data with theory/correlations, it is worthwhile to present the 
behaviour of stress and pressure in the packed bed so that hysteresis phenomenon can be 
understood properly. 

Equations (14) and (21) describe the stress distribution in both cartesian and radial 
regions respectively. Figure 6 shows the stress distribution as a function of distance in a 
packed bed above the cavity region. Stresses have been plotted from the top of the bed to 
the cavity roof for both the decreasing and increasing velocity at 40 m/s. In both the 
cases, the normal stress increases at a decreasing rate in the constant velocity region (i.e. 
z = 0.0m to z = 0.84m). Beyond this i.e. in the radial region, the stress keeps on 
increasing. The increment in radial stress continues until it reaches a maximum value, 
few centimeters away from the cavity roof. After reaching this maximum value, stress 
starts decreasing rapidly as it approaches the cavity roof. On closely scrutinising the 
equation (20) it can be deduced that pressure gradient term is responsible for this 
behaviour of the stress in the radial region. This can also be seen from Figure 7 in which 
pressure gradient has been plotted against the distance from the bed surface for both 
decreasing and increasing cases. Pressure gradient, given by Ergun equation (7), is 
constant in the cartesian region, however, pressure gradient in the radial region is a 
function of the distance from the center of the cavity and increases asymptotically close 
to the cavity roof. From Figure 7, it is clear that the value of pressure gradient is two 
orders of magnitude greater (close to the cavity roof) than the value of pressure gradient 
in the constant velocity region. 

The very high pressure gradient close to the cavity is responsible for the decrease in 
radial stress near the cavity roof. In the increasing velocity, the normal stress is always 
greater than the decreasing velocity as pressure drop is always higher in former case 7. 
This is one of the reasons that cavity size is less in the increasing velocity. 

In order to verify the proposed theory, experimental and published data have been 
compared with the theoretical predictions and results are presented below. 
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A comparison between the measured (published, Apte et aL (1990)) static pressure with 
the present theory is shown in Figure 8. As such the agreement between the measured and 
theoretical data is good except near the raceway roof which is due to the severe 
fluctuations in pressure which could lead the error in measurement up to 30%. Also a 
small difference in the location of measuring probe could give a very different result. 
Besides the above mentioned factors, there is also the discrepancies in the measured static 
pressure values as reported in the paper. For example, table in the paper shows raceway 
size 0.041m while figure shows 0.035m. Similarly, the bed height is reported 0.5m while 
figure shows 0.55m. However, for comparison we have taken operating data from the 
table as mentioned in the paper. Published pressure gradient values, Apte et al. (1990), 
(derived from the pressure measurement curve) also show a good agreement with the 
present theory which has not been reported here. 

From equation (28), it is clear that it can be solved for the cavity radius R as other 
parameters are known. Prediction of the cavity size in the increasing velocity is given in 
figure 9. As seen from this figure, when one starts increasing the velocity, cavity does not 
form until a critical velocity is reached. Initially, there is no cavity formation since the 
pressure force exerted by the gas is unable to overcome the frictional force and bed 
weight. Later when the velocity (or pressure force) is sufficiently high to overcome these 
forces, void or cavity starts forming. From this point onwards as one increase the 
velocity, the cavity size keeps on increasing. In another way, one can say that when the 
bed is in static condition (no gas flow), friction is acting in upward direction (due to the 
downward movement of the particles during bed filling). As soon as the gas is introduced 
friction try to resist the force exerted by the gas. As the pressure force keeps on 
increasing with the increase in velocity, frictional force starts acting in reverse direction 
and becomes fully mobilized when the formation of cavity occurs. Figure 9 also shows a 
comparison between experimental (Sarkar et al., 2003) and theoretical values of cavity 
size. Within the experimental errors, very good agreement exists. 

Figure 10 compares the theoretical and experimental cavity hysteresis. The constant 
region of cavity size in theoretical prediction is plotted based on the following arguments. 
Theoretically, it was found that the normal stress (and thus the frictional force) was 
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higher in the velocity increasing case than in the velocity decreasing case at a particular 
gas velocity (see Figure 6). The frictional force at the maximum gas velocity in the 
velocity increasing case is known (from an almost similar equation as in decreasing 
velocity (26)) and thus the maximum cavity size. In the velocity decreasing case, this 
maximum cavity penetration was taken as constant for each gas velocity until the 
frictional forces in the decreasing case attained a value equal or lower than the frictional 
force corresponding to the maximum gas velocity in the velocity increasing case. In other 
words, it is assumed that the frictional forces in the velocity decreasing case are fully 
mobilized when they become almost equal to the frictional force corresponding to the 
maximum gas velocity in the velocity increasing case. Therefore, one keeps the cavity 
size constant in the decreasing velocity until it obtains a value lower than that obtained at 
the maximum blast velocity in the velocity increasing case. Based on the above 
arguments, one can present the results in the region where there is not much change in the 
cavity penetration with the blast velocity as shown in Figure 9. A reasonable agreement 
can be seen at decreasing gas velocity between the two values. Slight disagreement in the 
value could be due to two reasons. Firstly, during experimental measurement of cavity, 
its size can vary by ± two particles diameter. The difference between the theoretical and 
experimental values is not more than two particle diameters as it can be seen from figures 
9 and 10. Secondly, the values of internal and wall angles of friction may change with 
change in pore pressure as reported by Terzaghi, K., Peck, R.B. & Mesri, G. 1996 Soil 
Mechanics In Engineering Practice. 3 rd Edition John Wiley New York, and Jackson and 
Judd (1981). A constant value of these angles at all gas flow rate has been considered in 
the present study. Nevertheless, it can be concluded that using the frictional force in the 
force balance, it is possible to predict the cavity penetration in the velocity increasing and 
decreasing cases with a fair degree of accuracy. 

Figure 11 shows a plot of cavity diameter vs. gas velocity without considering the 
frictional forces along with a theoretical hysteresis curve. It is obvious that a force 
balance approach based on gas momentum and bed weight terms can not give the correct 
results. Certainly, frictional forces play an important role in describing the packed bed 
behaviour. Moreover, these two forces can not explain the hysteresis phenomena as they 
will give only one set of data in both the cases i.e. increasing and decreasing gas velocity. 
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Discussion, Novelty and Inventive Steps: 

Two important observations were made during this study. Firstly, frictional forces play an 
important role in the overall force balance to describe the hysteresis phenomena as shown 
in Figure 1 1 . Only after considering these frictional forces one would be able to describe 
the behaviour of a packed or fluidized or spouted bed properly. Tsinontides, S.C. and 
Jackson, R. 1993 The mechanics of gas fluidized beds with an interval of stable 
fluidization. J. Fluid Mech. 255, 231 -21 A, had neglected the wall-particle frictional force. 
However, they had recognized the importance of it but were unable to take into account 
of it in explaining the results. Indeed, it has been found from the theory, presented here, 
that as the angle of wall-particle is decreased, hysteresis in the bed reduces (and thus the 
cavity size increases) which is in agreement qualitatively with the experiments reported 
by Sarkar et al. (2003). In fact, if friction were removed completely, we would argue 
there should be no hysteresis at all. Obviously, friction has a pronounced effect on 
hysteresis and can not be neglected. Secondly, there have been wide variations in the 
literature about handling the direction of stresses as explained earlier in the prior art 
section. This we have clearly demonstrated by theory, presented here, and experiments 
done by Sarkar et al. (2003), that frictional forces will change the direction depending 
upon the upward or downward movement of solids. From this theory, it is clear that while 
increasing blast velocity, the frictional force will act downward. However, Apte et al. 
(1990) have considered it in the upward direction. In case of Doyle III et al. (1986), since 
there is a downward movement of the solids, the shear stress will act in the upward 
direction. They have considered it in the downward direction. Pressure force always acts 
in the upward direction and bed weight always acts in the downward direction. Indeed, 
Chong, Y.C., Teo, C.S. & Leung, L.S. 1985 Encyclopedia of Fluid Mechanics, 
Cheremisinoff, N.P. (Ed) 4, 1127-1144, had also mentioned about it in explaining the 
hysteresis in fluidized bed. Tsinontides & Jackson (1993) tried to explain fluidization and 
bubbling regime hysteresis by introducing a complicated assumption of compressive and 
tensile stresses which they had related to void fraction. As such there is no significant 
variation in the void fraction in a packed bed, therefore, the above concept may not be 
applicable. 
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Two new correlations have been developed to predict the cavity size in increasing and 
decreasing gas velocity in the stationary/moving bed. These correlations were developed 
based on a systematic experimental study and then applying dimension analysis taking 
care of frictional properties of the particulate. No one has done this study before and 
developed the correlations in a systematic way as discussed in prior art section. Similarly, 
one-dimensional analytical mathematical model has been developed based on force 
balance approach proposed by us as new theory and discussed in the prior art section. 
Again, no investigator has developed such a model. Also no mathematical model is 
available which can describe the hysteresis phenomena in packed/spouted/fluidized beds 
except the one which has been developed here. Both the correlations and mathematical 
model have tremendous industrial application potential some of them have been 
discussed under examples section below. 

Brief Description of the Accompanying Drawings: 

In the drawings accompanying the specification, 

Figure 1 illustrates the experimental cavity hysteresis with the packed bed. 

Figure 2 illustrates the packed bed showing the essential regions used for modeling. 

Figure 3 illustrates the forces acting on an element in the cartesian region. 

Figure 4 illustrates the forces acting on an element in the radial region. 

Figure 5 illustrates the schematic diagram of the experimental setup. 

Figure 6 illustrates the hysteresis curve of normal stress at velocity = 40 m/s. 

Figure 7 illustrates the variation of pressure gradient with distance from the bed surface. 

Figure 8 illustrates a comparison of static pressure with experimental data. 

Figure 9 illustrates a comparison between theoretical and experimental cavity size for 

increasing velocity. 

Figure 10 illustrates a comparison between theoretical and experimental cavity 
hysteresis. 

Figure 11 illustrates a comparison between theoretical cavity size considering and not 
considering the frictional forces. 

Figure 12 illustrates a comparison of correlation raceway with published Flint and 
Burgess (1992) data of 3mm polystyrene. 
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Figure 13 illustrates a comparison of correlation raceway with published Flint and 
Burgess (1992) data of 0.725mm ballotini glass. 

Figure 14 illustrates a comparison of model's prediction with experimental (Born, 1991) 
values of cavity size. 

Figure 15 illustrates a comparison of experimental (Sastry, 2000) and theoretical values 
of cavity size. 

Figure 16 illustrates a comparison of raceway size between experimental and correlation 
in both increasing and decreasing velocity. 

Figure 17 illustrates a comparison of blast furnace (Hatano et al., 1977) and experimental 
data for both increasing and decreasing velocity conditions. 

Figure 18 illustrates a comparison of correlation raceway size with published (Wgastaff, 
1957) blast furnace data. 

Figure 19 illustrates a comparison of correlation raceway size with published blast 
furnace data of Nishi et al., 1982. 

Figure 20 illustrates a comparison of correlation raceway size with published blast 
furnace data of Poveromo et al., 1975. 

Figure 21 illustrates a flow chart for determination of cavity/raceway size in packed bed 
like ironmaking & lead blast furnaces, corex, cupola, etc. for decreasing gas velocity 
based on mathematical model. 

Figure 22 illustrates a flow chart for determination of cavity/raceway size in packed bed 
like ironmaking & lead blast furnaces, corex, cupola, etc. for decreasing gas velocity 
based on decreasing correlation. 

Figure 23 illustrates a flow chart for determination of maximum velocity/cavity size in a 
spouted bed above which spout will form / or condition of instability in packed bed based 
on mathematical model. 

Figure 24 illustrates the flow chart for determining cavity/raceway size in packed bed 
like ironmaking & lead blast furnaces, corex, cupola, etc. 

EXAMPLES: 
Examples 1-5: 

It was discussed in the beginning that many correlations have been proposed to predict 
the cavity size but they are not in agreement with each other. Now, it is the time to 
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validate the proposed correlations and mathematical model to see whether they can 
represent the experimental data of other researchers'. 

Figure 12 shows a comparison of raceway diameter obtained using the correlation and 
published experimental values for a 2D cold model (Flint and Burgess (1992)). The 
experimental values of the raceway diameter have been obtained for polystyrene beads of 
diameter 3mm, bed height from the tuyere level 800mm, and tuyere opening 5mm. Angle 
between the wall and particle was taken 18 (F. Bom, B.E. (Hons) Thesis, University of 
Queensland, Australia, 1991). Other values are given in Flint and Burgess (1992). 
Average raceway diameter was used in plotting the value as data were available for 
raceway penetration and raceway height. There is good agreement between the 
experimental values of average raceway diameter and that obtained using the correlation 
with the maximum error equal to the two particles diameter except at the maximum blast 
velocity which is close to the fluidisation limit and one can not expect a good agreement 
for those values. It should be noted that correlation gives the diameter of the raceway 
which could be little different from the raceway penetration or average raceway diameter 
as reported by many investigators. In the same figure, prediction by mathematical model 
has also been plotted and one can see an excellent agreement between the experimental 
and theoretical results. 

Another published experimental result (Flint and Burgess (1992)) along with correlation 
data for the glass bead of diameter 0.725mm, bed height from the tuyere level 800mm 
and tuyere opening 5mm is given in Figure 13. The results are for 2D cold model in 
increasing gas velocity condition. Average diameter of the raceway has been calculated 
from the published value and has been plotted in this figure. 

Wall-particle angle was taken 12.4 (Apte et al. (1990)). There is good agreement between 
the experimental values of average raceway diameter and that obtained using the 
correlation. 

Born (1991) has used polystyrene beads for 2D experiments using apparatus number 1 
(see Table 2). Comparison of model's predictions with Bom (1991) experimental data, in 
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increasing velocity, is shown in Figure 14. Because these data were used to develop the 
correlations therefore, correlations results are not compared. A good agreement is 
apparent between the two. 

G.S.S.R.K. Sastry, M.Sc. (Engg) Thesis, Indian Institute of Science, Bangalore, Sept. 
2000, has used quartz as particulate material in his increasing velocity experiments using 
2D apparatus number 4 (see Table 2). The properties of these materials are given in Table 
3 along with other parameters. A comparison is shown in Figure 15. Correlation's results 
are not shown in this figure because these data were used to develop the correlation. 
Again a good agreement is apparent between the two. 

Our experimental data has already been compared with mathematical model in Figure 10 
above. 

A comparison between the experimental and predicted (using correlation equation (35)) 
raceway size for decreasing gas velocity is shown in Figure 16. The plot is for plastic 
beads of diameter 2.1mm. Bed height was 600mm from tuyere level and tuyere opening 
was 5.5mm. Apparatus number 3 (see Table 2) was used during the experiments. An 
excellent agreement between them is apparent. 

It should be mentioned here that plastic beads data were not used to develop the present 
correlations. The linear decrease in the raceway penetration with blast velocity is 
predicted well by the correlation. Similarly, one can see a good agreement between the 
two values in increasing gas velocity as shown in the same figure. Mathematical model's 
results are also shown in the same Figure 16 for increasing velocity. A good agreement is 
evident. 

Examples 6-8: 

It was mentioned in the prior art section that raceway size obtained in decreasing gas 
velocity is more relevant to operating blast furnaces than increasing gas velocity. It is 
because large amount of coke is consumed near the raceway during combustion and in 
reducing the pre. This coke is replenished from the top of the raceway. Also 
intermittently iron and slag is tapped from the bottom due to which coke descends. It has 
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also been found (MacDonald & Bridgewater, 1993) that the decreasing gas velocity 
condition is applicable to the case of a moving bed as in the case of blast furnace. It was 
observed that the horizontal injection into a moving bed gives effects similar to those 
encountered with vertical injection into a moving bed. So the decreasing correlation 
results can be applied to the moving bed irrespective of whether there is horizontal or 
vertical injection of the gas. 

All the previous correlations, which have been given for the raceway penetration till now, 
are mainly for the increasing velocity. There is a doubt of their applicability to the blast 
furnaces. Now, it is the time to verify two points: 

1 . Whether the decreasing gas velocity is relevant to blast furnace or increasing, and 

2. Whether the developed correlation, based on cold model results, can represent the 
commercial blast furnace. 

Figure 17 shows the raceway factor values plotted as a function of penetration factor 
obtained using the correlation in the both velocity increasing and decreasing cases along 
with the blast furnace data (HatanO et al., 1977). For the sake of comparison, data 
obtained from mathematical model have also been included in the same figure. In this 
figure data obtained from the correlation and model are based upon the cold model 
experimental data for the case of apparatus 1, tuyere diameter 6mm, bed height lm and 
polyethylene beads of equivalent diameter 4.1 mm. There is an excellent agreement 
between the raceway factor values obtained in the velocity decreasing case with the blast 
furnace data when plotted as a function of penetration factor. It confirms both the points 
mentioned above that raceway size obtained in decreasing velocity is more relevant to 
commercial blast furnace and the correlations/mathematical model developed here 
reasonably predict the raceway size. 

In Figure 17, it was difficult to compare the raceway size against gas velocity using 
actual blast furnace data due to unavailability of many data. In fact, in most of the 
published work on commercial blast furnace, many data are missing. However, we have 
managed to extract most of the data from these papers. Some of the values have been 
assumed in a reasonable way which are described during the discussion of a particular 
figure. 
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Wagstaff et ah (1957) reported the data of commercial blast furnaces almost half a 
century ago. That time blast furnace technology was not so advanced. We were able to 
extract most of the data which are required by the correlation and model to predict the 
raceway size except the height of the burden, coke size, apparent density of solid and 
hearth radius (as W in the correlation). After going through few text books (The Iron and 
Steel Institute of Japan, Blast Furnace Phenomena and Modelling, Elsevier Applied Sci., 
London, (1987) and A.K. Biswas: Principles of Blast Furnace Ironmaking, SBA 
publications, Calcutta, 1984) coke size was assumed 40mm and apparent density of coke 
was assumed 900 kg/m 3 . These values were kept constant in other papers also (if 
applicable). Hearth diameter, especially for the old furnaces of 1950, was assumed 7m. 
Burden height was calculated, for all authors, as effective burden height using formula 
suggested by Sastry et al. (2003). They have shown, based on stresses at the bottom of a 
2D apparatus and using modified Janssen equation for two-dimension case, that pressure 
becomes almost constant at the bottom of the apparatus after certain burden height. Using 
their formula and assuming 15m burden height, it was found that pressure at the bottom 
becomes constant after 5m of burden height. Therefore, this height (5m) was taken as 
effective burden height for all commercial furnaces. It was also found that if burden 
height is taken to 20m then there is hardly any change in the effective bed height. 

Figure 18 shows a raceway size comparison between the plant data (Wagstaff et al., 
(1957)) and correlation with gas velocity. Correlation data have been plotted for 
decreasing velocity. Similarly, decreasing data obtained from the model has also been 
plotted in the same figure. Errors bar are also shown in the plant data. It is pleasing to see 
an excellent agreement between them. Because our correlations and mathematical model 
are based on two-dimensional model, therefore, tuyere diameter area was converted in to 
equivalent 2D tuyere area and then Dt, tuyere opening, was calculated. Diameter of the 
furnace tuyere was taken as the thickness of the apparatus for slot tuyere. 

Figure 19 shows another comparison between the correlation and Japanese blast furnaces 
(T. Nishi, H. Haraguchi, Y. Miura, S. Sakurai, K. Ono and H. Kanoshima, ISIJ, 1982, 
vol. 22, pp. 287-296). In this paper all the data were available except apparent density of 
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coke which was taken 900 kg/m 3 as described before. Again a good agreement exists 
between the two. The difference between the two values is mostly within the limit of ± 
two to four particles diameter. For comparison purpose, increasing velocity data, has also 
been plotted in the same figure. It is obvious that decreasing velocity data are relevant to 
blast furnaces and are well represented by decreasing cold model correlation. 

Another comparison of correlation with operating blast furnace data (J. J. Poveromo, 
W.D. Nothstein and J. Szekely: Ironmaking Proc, 1975, vol. 34, pp. 383-401) is shown 
in Figure 20. In this figure also, almost all data were given in the literature except coke 
size and its apparent density which were taken as 40mm and 900 kg/m 3 respectively. 
Again a good agreement exists between the two. Figure also shows the errors bar in the 
plant data. 

Example 9: 

The model developed here has provided a basic frame work to describe the complex 
phenomena of hysteresis in packed, fluidized and spouted beds including the stresses 
(between the particles and wall and particles) in a force balance which include gas drag 
and particles weight. At this point, it is important to make some comments on the nature 
of the equation (21). Stress can be estimated using equation (21). From this equation it 
can be seen that a r is strongly dependent on the pressure drop in the bed. Under fluidized 
bed condition, the bed weight is equal to pressure drop and thus a r would be zero. If 
pressure drop is greater than bed weight then a r may become negative. However, in the 
packed bed, particles are in contact with each other and with the container wall therefore, 
a r may not achieve a negative or zero value unless the bed approaches fluidized 
condition. This is an important conclusion as Apte et al. (1990) assumed that a r could 
achieve a negative value above the cavity roof. Tsinontides & Jackson (1993) has also 
ruled out that a r could achieve a negative value. Obviously, Apte et al. assumption was 
incorrect in explaining experimental hysteresis. Using equation (21), the velocity at 
which a bed may become unstable/fluidized can be found provided all the properties of 
the particulate material and gas are known. From Figure 6 it is clear that there is a 
maximum value of radial stress exist above the cavity roof for a particular gas velocity. 
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This maximum value of stress decreases with increase of gas velocity. This indicates that 
to make the bed unstable/fluidized, one has to overcome this maximum stress by 
increasing the gas velocity. Therefore, it is possible that at a particular gas velocity the 
system can overcome this maximum value of stress and will become unstable/fluidized. 
Indeed, during our experiments it was found that the bed becomes unstable when we 
approach a velocity close to 142 ± 5 m/s. Using equation (21), the velocity for unstable 
bed was found 140 m/s which is an excellent agreement. 

Conclusions: 

Two raceway size correlations have been developed one each for increasing and 
decreasing velocity under the cold model conditions. Frictional properties of the material 
have also been included in these correlations. Raceway size obtained from the 
correlations and other data such as published cold & hot model, plant and experimental 
data match very well. It has been shown that decreasing conditions prevails in the 
operating blast furnace and therefore, decreasing correlation can be used to predict the 
raceway size. Both the correlations are able to predict the raceway hysteresis in cold 
model. It has been found that the frictional forces (and thus the frictional properties) have 
pronounced effect on the prediction of cavity size. In fact, the inclusion of frictional 
forces gives a universal form to the force balance approach to predict the cavity size. This 
is evident by comparing the theoretical data with published, experimental and plant data. 
An excellent agreement has been found between theory and experiments, not only with 
our experiments but also with other researchers experiments under various conditions. 
With the help of mathematical model, the maximum operating velocity of any packed bed 
can be found, above which the bed may become unstable and thus its operation. 

The main advantages of the present invention are: 

1. It can explain the hysteresis phenomena correctly and shows the importance of 
frictional forces in packed/fluidized/spouted bed systems. 

2. It can predict the cavity size in these systems so that their performance can be 
improved considerable in terms of heat, mass and momentum transfer. 

3. It gives two simple working correlations, beside a mathematical model, to predict 
the cavity size one each in increasing and decreasing velocity respectively. 
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4. It shows cavity size pertaining to decreasing gas velocity is relevant to operating 
blast furnaces. 

5. Mathematical model can also give the maximum operating gas velocity for any 
packed bed above which it may become unstable. 

6. Both model and correlations have been tested under wide variety of conditions 
(see examples 1-9) and they give very good results. Therefore, they can be used 
directly in the industries. 

7. No other models and correlations have the above mentioned features until now. 
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